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Scanning tunneling spectroscopy measurements on the high temperature superconductor 
Bi2Sr2CaCu208+i have reported an enhanced spectral gap in the neighborhood of O dopant atoms. 
We calculate, within density functional theory (DFT), the change in electronic structure due to 
such a dopant. We then construct and discuss the validity of several tight binding (TB) fits to the 
DFT bands with and without an O dopant. With the doping-modulated TB parameters, we finally 
evaluate the spin susceptibility and pairing interaction within spin fluctuation theory. The d-wave 
pairing eigenvalues are enhanced above the pure system without O dopant, supporting the picture 
of enhanced local pairing around such a defect. 



PACS numbers: 

I. INTRODUCTION 

With the observation in scanning tunnehng spec- 
troscopy (STS) measurements on high-Tc Bi-based 
cuprate superconductor^^ that the oxygen dopant posi- 
tion and the size of the superconducting gap are corre- 
lated, it has become evident that dopant atoms may in- 
fluence superconductivity beyond their roles as sources of 
mobile charge and scattering centers. Subsequently, Nun- 
ner et aL^^l showed that many STS observations can be ex- 
plained if one assumes that dopants locally enhance the 
pairing interaction. It is clearly desirable to further jus- 
tify the proposal of Ref. [2]by identifying the microscopic 
mechanism responsible for the local pairing interaction 
enhancement. Such a step would fulfill a longstanding 
goal of allowing the systematic study of the pairing in- 
teraction itself by the measurement and modelling of the 
effects which modulate it. 

If we accept the dominance of a magnetically driven 
pairing mechanism in the cuprates, the pairing enhance- 
ment is a result of the local increase in the spin fluctu- 
ation exchange interaction in the vicinity of a dopant, 
which can occur due to the dopant-induced local struc- 
tural modifications. Whether the structural and the cor- 
responding local electronic structure modifications can 
indeed enhance the local superexchange coupling has 
been a subject of several recent studies'^ The prob- 
lem is typically treated by calculating the local exchange 
coupling constants of the t-J model, the large-[/ limit 
of the Hubbard model, with the assumption that the lo- 
cal parameters of the latter (the on-site energy fi and the 
hopping integrals t) are modified due to the presence of a 
dopant. In Ref. by means of the Rayleigh-Schrodinger 
perturbation expansion it was shown for the three-band 
Hubbard model that the exchange coupling is enhanced 
only in a certain region of the Hubbard model parame- 



ters phase diagram. Unfortunately, recent numerical esti- 
mates of the parameter values of the Hubbard model with 
an impurity based on electrostatic calculations^ place 
the dopant-induced variation in the region of the phase 
diagram where exchange gets suppressed. Though dis- 
couraging at first glance, it is clear that this general 
approach contains many oversimplifications of the true 
electronic structure. In addition, the full pairing inter- 
action, while related to the exchange coupling J, can be 
influenced also by dynamical spin fluctuation processes.'^ 
We have therefore been stimulated to investigate further 
the dopant-induced effects on the local pairing with the 
goal of improving the model by a reflnement of the ap- 
proximations. For instance, besides considering the vari- 
ation of the atomic on-site energies (which was assumed 
in Ref . El and E] to be the only effect due to a dopant) one 
can allow for the variation of hopping integrals near the 
dopant as well. One can furthermore go beyond the elec- 
trostatic considerations in calculating the inhomogeneous 
Hubbard model parameters in order to place more accu- 
rately the exchange coupling variations in the parameter 
phase diagram by employing ab initio calculations. Fi- 
nally, one can directly calculate the dynamical pairing 
interaction within weak coupling spin fluctuation theory. 
These are the aims of the present work. 

In this paper, we perform DFT calculations for 
Bi2Sr2CaCu208+5 and explore the possibility of extract- 
ing out of the electronic structure reliable effective Cu-Cu 
hopping integrals t and on-site energies /i as a function of 
oxygen doping. In general, DFT calculations of the elec- 
tronic structure near the Fermi surface (FS) of cuprate 
materials are not considered reliable due to strong corre- 
lation effects, but it has been argued that changes in elec- 
tronic structure induced by high energy impurity states 
are much less sensitive. These effective parameters are 
further considered in the Hubbard model calculations of 
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the spin susceptibiUty and superconducting gap func- 
tion, and key changes in these functions are observed 
which support the importance of the effect of the oxygen 
dopants on the pair correlations in these materials. 

The evaluation of the hopping integrals and on-site 
energies of the Hubbard model is performed by map- 
ping the eigenvalues of the corresponding non-interacting 
tight-binding Hamiltonian to the valence bands of 
Bi2Sr2CaCu208+5 as a function of oxygen doping. Al- 
though Bi2Sr2CaCu2 08+5 is more accurately described 
by the three-band model^, here we consider a single-band 
TB Hamiltonian, with the single active orbital being the 
Cu 3dx2-y2 orbital at the Fermi level as the simplest 
case for the proposed approach. In fact, recent studies 
showed that the one-band Hubbard model suffices to de- 
scribe the irnportant features of the low-energy physics in 
the cupratea^. We find that the accuracy of the dopant- 
induced TB Hamiltonian parametrization suffers from 
two sources of uncertainty, which are (i) the presence of 
effective far Cu neighbor interactions in the homogeneous 
model (undoped case) allowing for many alternative pa- 
rameter sets and (ii) the need to reduce the enormous 
number of adjustable parameters of the dopant-induced 
TB model in order to perform the numerical optimiza- 
tion. In view of these considerations we present two al- 
ternative TB models and discuss their validity in terms 
of physical arguments. 



II. ELECTRONIC STRUCTURE 
CALCULATIONS 

The electronic structure calculations for the par- 
ent compound Bi2Sr2CaCu20g were performed with 
the reference crystal structure reported in Ref. [S] 
Bi2Sr2CaCu208 crystallizes in the space group lA/mmm, 
with a unit cell which we consider in the tetragonal sym- 
metry (see Ref. ^ for the discussion on the structural 
supermodulation in Bi2Sr2CaCu208), consisting of two 
identical slabs of atoms, one shifted with respect to the 
other by a vector (a/2, a/2, c/2), where a and c are the 
lattice parameters of a tetragonal unit cell. Oxygen dop- 
ing of this system was modeled by introducing one extra 
interstitial oxygen atom into a surface supercell consist- 
ing of eight primitive unit cells in the xy plane. We 
consider the surface supercell in order to reproduce the 
conditions of a typical STM experiment on an 0-doped 
Bi2Sr2CaCu208+5 surface. The electronic structure cal- 
culations in the present work were performed with the 
linearized-augmented-plane-wave (LAPW) basis, as im- 
plemented in Wien21i[12J. The exchange and correlation ef- 
fects were treated within the generalized gradient approx- 
imation (GGA) as implemented by Perdew, Burke and 
Ernzerhof23. Additional details of the electronic struc- 
ture calculations are given in Appendix A. 

In Fig. [l] (a) we present the electronic bandstruc- 
ture of Bi2Sr2CaCu208 in the energy window between 
[—8 eV, 3 eV]. The two Cu 3dx2-y2 bands of the parent 



compound are marked by circles with sizes proportional 
to the 3dx2-y2 weight. They are rather dispersive, which 
is a typical feature of the cuprate familjff^l xhe two Cu 
bands crossing the Fermi level are the anti-bonding bands 
in the chemical bonding between Cu atoms and the O 
atoms in the superconducting layer. The bonding bands 
lie in the region between [—8 eV, —5 eV]. In the follow- 
ing we concentrate on the energy window [—2.5 eV, 3 eV] 
near the Fermi level Ep as indicated in Fig.[l](b). At the 
X = (7r,0, 0) point the Cu 3dx2-y2 bands show a strong 
overlap with the Bi-0 bands near the Fermi level, and 
at the r = (0, 0, 0) point there is strong hybridization 
between the Cu 3dx2-y2 and some of the lower- lying O 
2p and Cu 3d bands. 

Fig. [T] (c) presents the bandstructure of the 0-doped 
Bi2Sr2CaCu208+5 supercell in the same energy range. 
The bandstructure of the supercell is plotted along the 
same high symmetry points in the Brillouin zone as those 
of the parent compound. However, since the doped sys- 
tems' unit cell is 8 times larger than the parent com- 
pound unit cell, its Brillouin zone shrinks and the bands 
fold on top of each other so that effectively there are 16 
Cu 3dx2_y2 bands crossing the Fermi level, where 16 is 
the number of Cu atoms in the supercell slab. 

Because the Brillouin zones of the parent compound 
and of the doped supercell are defined differently, it is 
not straightforward to compare their bandstructures. In 
order to make such a comparison, we recalculate the elec- 
tronic structure of the parent compound in the folded 
Brillouin zone. Once this is done, small but traceable 
changes in the shape of the Cu 3d.j;2_y2 bands as a func- 
tion of oxygen doping become apparent as will be shown 
in the next sections. We will dedicate the rest of the 
paper to the quantitative evaluation of these changes in 
terms of hopping integrals of a single-band tight-binding 
Hamiltonian as well as calculations of the spin suscepti- 
bility and superconducting pair function. 



III. SINGLE-BAND TIGHT-BINDING MODEL 
A. Parent compound 

A standard procedure to analyze the low-energy band- 
structure features of solid st ate s ystems is the tight- 
binding (TB) parametrizatiorP^li^. Since the symme- 
try of the parent compound Bi2Sr2CaCu208 unit cell is 
lA/mmm and the number of distinct hopping integrals 
(adjustable parameters) is small, the TB parametrization 
in this case is straightforw ard a nd has been already per- 
formed in previous studieJi^'^. The complexity of the 
problem is dramatically increased when trying to obtain 
the TB parameters for the doped supercell bands shown 
in Fig. [I] (c): first, because the number of bands to be 
mapped increases by a factor of 8, which means that, 
in order to optimize the TB parameters, a global mini- 
mum of a complex mathematical function expressed by 
a 16 X 16 matrix needs to be found; second, because the 
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FIG. 1: (Color Online) (a) and (b) DFT electronic bandstructure of Bi2Sr2CaCu208 shown for different energy ranges, (c) 
Bandstructure of the O-doped Bi2Sr2CaCu2 08+i supercell in the same energy window as (b). The weight of the Cu 3d^2_y2 
character in the eigenvalues is proportional to the circles' size. 



number of distinct hopping integrals is expected to rise 
considerably as we increase the size of the unit cell and 
lower its symmetry. Our strategy to overcome these com- 
plications is to use the hopping integrals obtained for the 
parent compound Cu Sd^^-y^ bands as starting values 
for parametrizing the doped supercell bands. 

An important point to consider is that the differences 
between the Cu idx^-y^ bands of the parent and doped 
compounds are small, which is reasonable since we do 
not expect the interstitial oxygen to have a drastic ef- 
fect on the orbital overlap of its neighboring Cu atoms. 
Therefore, the TB parameters for the Cu 3dx2-y2 bands 
in the parent and doped compounds must be obtained 
with a fitting error smaller than the dopant-induced dif- 
ferences in the bandstructure, i.e., the quality of both 
TB mappings must be particularly high. For this rea- 
son we consider the antibonding Cu id^i^yi bands over 
the entire energy range over which they disperse. We 
find that, in order to accurately reproduce the DFT Cu 
3dx2_y2 bands of the parent compound [Fig. fll (a)], up 
to 13 Cu-Cu neighbors have to be included in the model 
TB Hamiltonian. We should note, however, that among 
these various hopping integrals, those representing hy- 
bridizations between the far neighbors are only effective 
parameters; we are forced to include them in the single- 
band model in order to describe those features of the Cu 
id.j.2^y2 bands that stem from the interaction of the Cu 
'idx2_y2 orbital with other Cu 3d orbitals and the O 2p 
orbitals. For example, near F = (0,0,0) the Cu 'idx2-y2 
bands are anomalously flat and show an Ck oc A:^ behav- 
ioi^, where k is the momentum in the Brillouin zone. In 
order to reproduce this feature, inclusion of higher har- 
monics are required in the model equations. 

The long range effective hoppings which arise from the 
mapping of the complex band structure over a ~ 5 eV 
range onto a single band model are found to be much 
smaller than the short range hoppings, e. g., ^loo and 
tiiQ. This strongly suggests that it will be difficult to 



identify a unique parameter set. Alternative single-band 
models can exist corresponding to different choices of ef- 
fective hopping paths. We therefore discuss two possi- 
ble sets of single-band TB Hamiltonian parameters (pre- 
sented in Tables [T| and [ll] of Appendix B) to give a sense 
of how robust the TB models can be. Details of the 
construction of these models are found in the Appendix. 
Both sets are indeed found to describe the Cu 3dx2_y2 
bands in the parent compound equally well. In the fol- 
lowing we will adopt the notation 

• TBI: 12 hopping parameters, 6 effective interlayer 
hoppings 

— TBlundopcci: Undoped parent compound. 

— TBlioc. doped: TB parameters are taken to 
vary in the vicinity of O dopant. 

• TB2: 13 hopping parameters, 4 effective interlayer 
hoppings 

— TB2undopcd: Undoped parent compound. 

— TB2hoin. doped: Fit is obtained with homoge- 
neous TB parameters over the entire supercell. 

A discussion of the tight-binding parameters describing 
the O-doped system will be the subject of the next Sub- 
section. 

In Fig.^(a) the TBI undoped Hamiltonian spectrum in 
the full Brillouin zone is compared to the DFT bands of 
the parent compound. In Fig. [2] (b) , the spectrum of the 
TBlundoped Hamiltonian is replotted in the folded Bril- 
louin zone and compared with the Cu 3dj:2-y2 bands of 
the doped supercell. The slight differences between the 
shape of the doped supercell DFT bands and of the TB 
bands of the parent compound are due to the presence of 
the interstitial oxygen, which displaces the neighboring 
Cu atoms and thus modifies their corresponding over- 
lap integrals. It is possible to quantify these effects by 
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FIG. 2: (Color Online) Plot of the DFT calculated Cu 3d^2_y2 
bands (dots) and the TBI Hamiltonian spectrum (lines) for 
Bi2Sr2CaCu2 08 and Bi2Sr2CaCu2 08+i : (a) comparison of 
the Bi2Sr2CaCu208 electronic structure to the TBlundopcd 
model; (b) comparison of the Bi2Sr2CaCu208+« electronic 
structure to the TBlundopcd model plotted in the folded Bril- 
louin zone; (c) comparison of the Bi2Sr2CaCu2 08+i elec- 
tronic structure to the TBlioc. doped model (see Fig.[4|. High 
symmetry points are given by {kx,ky) only, fez — 0; thus 
(0,0) = r, (^,0) =X, (tt,^) =T. 
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FIG. 3: (Color Online) Plot of DFT calculated Cu M^2_y2 
bands (dots) and the TB2 Hamiltonian spectrum (line) for 
Bi2Sr2CaCu2 08 and Bi2Sr2CaCu2 08+i : (a) comparison of 
the Bi2Sr2CaCu208 electronic structure to the TB2undopcd 
model; (b) comparison of the Bi2Sr2CaCu2 08+« electronic 
structure to the TB2undopcd model plotted in the folded Bril- 
louin zone; (c) comparison of the Bi2Sr2CaCu2 08+5 electronic 
structure to the homogeneous TB2hom. doped model. ( Table [ll| 
second row). 



fine-tuning the parameters of the TB Hamiltonian such 
that the supercell bands are reproduced, as we wiU show 
below. 

In Fig. |3] (a), we now show the corresponding 
TB2undopod Hamiltonian spectrum in the full Brillouin 
zone in comparison to the DFT bands of the parent com- 
pound. In Fig.[3](b), the same TB bands are plotted in 
the folded Brillouin zone, together with the DFT bands 
of the doped supercell. Comparing the TBlundopcd and 
TB2un(iopod model parameters for the parent compound 
(Tables g and |ll] of Appendix B), we observe that the 
TB2undopcd model has several features that can be con- 
sidered an advantage in terms of the physics that the 
model implies; this model includes only four hopping in- 
tegrals between the Cu02 layers, tooz, tuz, i2iz and ^332, 



whose relevance can be justified either by the close prox- 
imity of the two Cu atoms (tooz) or by the presence of 
a Ca atom along the Cu-Cu connection mediating elec- 
tron hopping (iii2, t2iz and tsa^)- In the TBlundopcd 
model, on the other hand, the mechanism of some of its 
six inter-layer interactions is not as clear. Furthermore, 
one would rather expect the contribution of interacting 
far Cu neighbors within a Cu02 layer to be more impor- 
tant as considered in the TB2undopcd model. 

In the next section, we will discuss the results of 
the derivation of the TB Hamiltonian for the 0-doped 
Bi2Sr2CaCu208+5 supercell, obtained when either the 
parameter set in TBlundopcd (Table|l]) or TB2undopcd (Ta- 
ble [h]) are used as initial values for mapping the supercell 
DFT bands. 
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B. O-doped supercell 

In order to describe the DFT Cu 3dj.2_y2 bands of the 
doped superceU [Fig. [T] (c)] within a TB model, one has 
to construct a Haniiltonian similar to Eq. ( |B1[ ), which 
will be now represented by a 16x16 matrix, according 
to the number of Cu atoms in the supercell (8 atoms 
in each Cu02 layer) . Since the presence of the intersti- 
tial oxygen introduces inhomogeneities in the system, the 
number of distinct model parameters for the supercell is 
not defined by simply the number of parameters of the 
corresponding parent unit cell (which is 12 and 13 for the 
models of Table |l] and Table [Til respectively, plus the on- 
site energy /i), but increases considerably. For instance, 
even by taking into account the mirror plane symmetry, 
there are still 238 parameters in the supercell TB Hamil- 
tonian based on the TBlundopcd ansatz. Technically, it is 
impossible to find a unique and unambiguous set of pa- 
rameters by performing an optimization of such a huge 
number of parameters, especially since our aim is to cap- 
ture the slight differences between the bands of the parent 
compound and the doped compound [see Fig. [2] (b) and 
Fig. 3](b)]. One way to proceed would be to approximate 
the hopping integrals that become distinct in the super- 
cell due to the inhomogeneity introduced by the dopant 
by their average values. In this "averaged" homogeneous 
TB model for the supercell, there would be as many pa- 
rameters as in the corresponding model for the parent 
compound, and their optimization would be simple; an 
example is given in Table [TT| (TB2hom. doped row). With 
such an approach, however, the most interesting physics 
concerning local effects due to the dopant is left out. 

The exact knowledge of how the Cu on-site energies 
and the most relevant Cu-Cu hopping integrals tioo and 
tiio are modified near the dopant is very important for 
understanding the dopant-induced effects on the local 
spin superexchange couplin^i^lli^, which is related to the 
size of the local superconducting gap in cuprateJ^. There- 
fore, in order to be able to study the local variations in 
the model parameters, we propose the following approxi- 
mate treatment of the problem. We assume that the on- 
site energies and hopping integrals most affected by the 
dopant are those that are nearest to the dopant, and we 
concentrate on the largest TB model parameters, such 
as fi, tiQQ and tuQ. Then, the supercell TB Hamilto- 
nian is optimized by adjusting the selected parameters, 
while for the rest, i.e., the effective far Cu neighbor in- 
teractions, their initial values are preserved. An illus- 
tration of the procedure in the case of the TBlundopcd 
model is shown in Fig. [4] and is worked out in detail 
in Appendix C. It turns out that this approach gives 
sensible results when the TB supercell model is built 
upon the TBlundopcd modcl (Table |l]) and gives coun- 
terintuitive results when the TB2undopcd model for the 
parent compound (Table |ll]) is used. In Fig. [2] (c) we 
show the good agreement between the energy spectrum 
of this TBlioc. doped Hamiltonian with the DFT bands of 
Bi2Sr2CaCu208+5. 
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FIG. 4: (Color Online) The schematic lattice of dopant- 
displaced Cu atoms in the oxygen-doped supercell of 
Bi2Sr2CaCu2 08+(5. The smallest circle represents the inter- 
stitial oxygen atom, and the larger circles stand for Cu atoms. 
Darker color is used for Cu atoms in the Cu02 layer closest to 
the interstitial oxygen. The Cu-Cu bonds that correspond to 
the hopping integrals of the tioo type are represented by solid 
lines, the hopping integrals of the two type by dashed lines, 
and the hopping integrals of the t200 type by dash-dotted 
lines. The numbers over the bonds stand for the optimized 
values of corresponding hopping integrals of the doped super- 
cell TB model based upon the nearest-neighbors parent com- 
pound model (see text). The optimized value of the on-site 
energy of this model for the Cu in the next-nearest Cu02 layer 
(the light Cu atom symbols) is = 0.4445 eV. The six opti- 
mized on-site energy values for Cu in the nearest Cu02 layer 
(the dark Cu atom symbols) are ^ = 0.5757 eV, 0.5057 eV, 
0.5341 eV, 0.5151 eV, 0.4930 eV, and 0.5186 eV. p = 0.5757 
eV corresponds to the Cu atom which is displaced most by 
the dopant and is marked with an arrow, m labels the mirror 
plane. 



To conclude this section, we present the homogeneous 
TB model for the doped supercell TB2iiom. doped based on 
the TB2undopod model (see Table [lT| and Fig. 31(c)). Even 
though the homogeneous model does not reflect the local 
dopant-induced effects, it can still be useful as it provides 
the information on how the model parameters change on 
average. For instance, it is interesting to observe that the 
average ratio ino/iioo increases in the doped supercell 
compared to that in the parent compound (from 0.2097 
to 0.2249); this might suggest a possible increase of the 
superconducting transition temperature upon doping, in 
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analogy with the observation within the cuprate family, 
that the materials characterized by a larger tug /<ioo ratio 
have higher transition temperatureJ^^. Unlike the situ- 
ation with the mapping approach previously discussed, 
which aimed at capturing local physics, the parameters 
of the homogeneous Hamiltonian demonstrate the same 
behavior upon doping regardless of the parent compound 
TB model (TBlundopod and TB2undopcd) chosen as a 
starting point for mapping the doped supercell electronic 
structure. 



IV. SPIN FLUCTUATION PAIRING 

A. Spin susceptibility 

In the following we calculate the magnetic spin suscep- 
tibilities for the TB models obtained previously, namely, 
(i) the TBlundopod (Table |l]), (ii) the TB2undopod (Ta- 
ble |ll] first row) , (iii) the inhomogeneous doped super- 
cell model TBlioc. doped (Fig. 4]) and (iv) the homoge- 
neous doped supercell model TB2hom. doped (Table [ll] sec- 
ond row). The spin susceptibility is derived within the 
Matsubara Green's functions formalisrrpS! from the non- 
interacting Green's functions. In a general formulation, 
the spin susceptibility is a function of four orbital indices, 
(Xs)st , which, in the considered case of a single orbital 
but multiple atoms in a unit cell, refer to the orbitals on 
different atoms. For the non-interacting case, the spin 
susceptibility (xs)s? is equivalent to the charge suscepti- 



bihty (Xe)^^ iXsfs't 



(Xc)f 



Xst : and is given bjA^ 



x°aq,w) 



1 



q)) - fiE^im 



a^(k)aP*(k)ag(k-fq)at*(k-t-q) 



(1) 



In this expression, indices s, p, q and t refer to the 
N Cu atoms in the unit cell and run from 1 to TV while 
indices /i and ly distinguish the N eigenvalues £^(k) of 
the diagonalized TB Hamiltonian. The matrix elements 

(k) are the components of the eigenvectors of the TB 
Hamiltonian. The integration over the Brillouin zone has 
been replaced by a sum over a sufficiently large num- 
ber TVk of fc-points. f{E) is the Fermi-Dirac distribution 
function. In the following, we will focus on the static 
non- interacting spin susceptibility xs(q)i 



xs(q) = l^x7s{(i,^ 



0), 



(2) 



sp 



and examine its behavior in the four cases of interest 
along the main symmetry directions in the Brillouin zone. 

The static spin susceptibilities xs(q) of the parent 
compound calculated with the TBlundopod model of Ta- 
ble |l] (bold black line) and with the TB2undoped model 
of Table [h] (bold dashed line) are plotted in Fig. [s] (a). 
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FIG. 5: The static spin susceptibility of (a) the two parent 
compound TB models TBlundopod and TB2undopod and the 
homogeneous doped supercell model TB2hom. doped, plotted in 
the full Brillouin zone, and (b) the two parent compound TB 
models TBlundopod aud TB2undopod and the inhomogeneous 
TBlioc. doped and homogeneous TB2hom. doped doped supercell 
models, plotted in the folded Brillouin zone. 



The two susceptibilities show similar features with dou- 
ble peaks along (0, 0, 0) - (tt, 0, 0) and (tt, tt, 0) - (0, 0, 0) 
directions and a broad plateau at (tTjTt, 0). These simi- 
larities can be understood in terms of the fact that most 
important parameters in the two TB models (iioo, ^iiOj 
etc.) have close values. In this respect, it is not surpris- 
ing that the spin susceptibility calculated with the av- 
eraged TB parameters of the homogeneous Hamiltonian 
TB2hom. doped for Bi2Sr2CaCu208+5 (thin dashed line in 
Fig. [5] (a)) qualitatively reproduces the same behavior as 
the TB models for the parent compound. 

We next calculate the spin susceptibility with 
the inhomogeneous TB Hamiltonian for the doped 
Bi2Sr2CaCu208+5 supercell TBlioc. doped- Of course 
in this case the spin susceptibility must be calculated 
with the full supercell (16 x 16) Hamiltonian matrix 
and accordingly is defined in the folded Brillouin zone. 
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Fig. [5 (b) shows the spin susceptibihty calculated with 
the inhomogeneous TBlioc. doped model (thin black line) 
and, for comparison, the two parent compound suscep- 
tibilities replotted in the folded Brillouin zone (as in 
Fig. [5] (a), bold black and bold dashed lines). Within 
the inhomogeneous model, a pronounced peak in the spin 
susceptibility evolves at (tt, 0, 0) upon doping whereas in 
both parent compound susceptibilities this region is fea- 
tured by a shallow minimum in between two asymmet- 
rical peaks located at some distance from (tt, 0,0). A 
peak in the spin susceptibility of a non-interacting sys- 
tem can transform into a divergence indicating magnetic 
instabilities and possible ordering, when the interparticlc 
interactions are switched on. (7r,0,0) corresponds to a 
commensurate antiferromagnetic striped order with pe- 
riod 2\/2a, with stripes along the (110) direction of the 
parent compound unit cell. 



B. Superconducting gap function 

We consider now the models TBlundopcd a-nd 
TBlioc. doped in order to analyze the superconducting 
properties of the undoped and doped Bi2Sr2CaCu208+5, 
respectively. We calculate the pairing vertex by assuming 
that superconductivity in the high-T^ cuprates is driven 
by the exchange of spin and charge fluctuation^. The 
many-body effects of the Coulomb interaction are here 
treated within the random phase approximation (RPA). 

In order to calculate the pairing vertex, the 
RPA charge and spin susceptibilities, x^^^i'i:^) a-^^d 
X^^^(q, w), are required. They can be obtained from 
the non-interacting susceptibility x(q, to) in the form of 
Dyson-type equations as 



by the projected interaction vertex 



( RPAnp? ^ p 
VAc J st As 



VAc Juv \ ^ JwzA,s 



(3) 



and 



uvwz 



For a single-band model, only the diagonal and 
matrices' components are non-zero: 



(5) 



where U is the strength of the on-site intra-band 
Coulomb repulsion between electrons. The singlet pair- 
ing vertex is then given by 



r??(k,k',c.) 



1 



(6) 



ic/^Xc^P^(k-k',..)C/^ + lc/^ 



tq 



ps 



The scattering of a Cooper pair from the state (k, — k) 
to the state (k', — k') on the Fermi surface is determined 



r(k,k',a.) = 5]at(-k)<(k)r??(k,k',u;) 

stpq 

X <;*(k')a^';(-k'). 



(7) 



where indices v and refer to the eigenvectors of the TB 
Hamiltonian with the corresponding energy eigenvalues 
close to the Fermi level. As the strength of the pairing 
interaction is defined by a frequency integral of the imag- 
inary part of r(k, k'jw) weighted by it is sufficient 
to consider the real part of r(k,k',a; = 0) according to 
the Kramers-Kronig relation: 



°°d^ ^"^[^^'^'^''")U Re[r(k,k-,a.^O)]. (8) 



If the superconducting gap is decomposed into an am- 
plitude A and a normalized gap function g(k), the latter 
can be evaluated from the following eigenvalue equation 



/dk.' 1 



Here, 



r^y"™(k, k') = -Re [r(k, k', 0) + r(k, -k', o)] (lo) 



is the symmetric part of the full interaction and 
VFik) - |VkS,(k)| 



(11) 



is the Fermi velocity at point k on the Fermi surface. 
The largest eigenvalue A of Eq. (|9| determines the super- 
conducting transition temperature and its corresponding 
eigenfunction ^(k) has the symmetry of the gap. 

We have solved the eigenvalue problem ^ for the un- 
doped and doped Bi2Sr2CaCu208+5 models in the folded 
Brillouin zone of the supercell. The folded Brillouin zone 
has been considered in both cases in order to ensure that 
the eigenvalue equations are constructed under the same 
conditions, which is important when the resulting pair- 
ing strengths are compared. The calculations have been 
performed for the temperature T = 0.01 eV and we con- 
sidered Coulomb repulsion U values that range from 1.00 
eV to 1.66 eV. Note that these values represent renor- 
malized values of the Hubbard U appropriate for RPA 
treatments and are smaller than bare J/'a^^. 

We find that the doped Bi2Sr2CaCu208+5 model is 
characterized by a larger value of the pairing strength 
A compared to that of the undoped model. The pairing 
strengths for the two models are presented in Fig.|6](a) as 
a function of U. Below U ~ 1.5 eV the two A values are 
almost equal, but at larger U values the pairing strength 
for the doped model grows faster and diverges at t/ = 
1.65 eV. 

Fig. [g] (b) displays the gap functions g{k) of the un- 
doped and doped models, corresponding to the leading 
eigenproblem solutions A of Fig. [6] (a) at U = 1.64 eV. 
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FIG. 6: (Color Online) (a) The pairing strength A for the un- 
doped (triangles) and doped (circles) Bi2Sr2CaCu2 08+a TB 
models as a function of Coulomb repulsion U. (b) The su- 
perconducting gap function g(k) on the fc-point mesh at the 
Fermi surface of the undoped Bi2Sr2CaCu2 08 unit cell for 
the undoped (triangles) and doped (circles) models. The red 
(blue) color represents positive (negative) g(k) values, and 
the intensity of the color is proportional to the absolute value 
of (?(k). Half the nesting vector quest is shown by an arrow. 



The result from Eq. ([9| g(k) is defined on the mesh 
of fc-points at the Fermi surface of the folded Brillouin 
zone. In Fig. [6] (b), the fc-point mesh was unfolded to 
the Brillouin zone of the undoped compound unit cell 
in order to allow a comparison of the Fermi surface be- 
havior for the two systems with experiment. One should 
note that in the case of the doped supercell such an un- 
folding is, strictly speaking, not allowed and results in 



a tearing of the Fermi surface. Yet, since the symme- 
try lowering effects caused by a dopant are small, the 
unfolding in this case is a reasonable approximation. In 
particular, the unfolded way of presenting ^(k) allows us 
to observe that the symmetry of the undoped model g(k) 
is d^i^yi and that upon doping it is roughly preserved, 
though slightly distorted. We also note the characteristic 
reduction of the norm of the nesting wave vector (/nest, 
associated with the occurrence of the spin density wave 
state, in the doped case. This behavior is in agreement 
with ARPES results for BizSrzCaCuaOg+i ISll. 

The gap equation calculations presented here 
show that the TB model derived for the 0-doped 
Bi2Sr2CaCu2084.5, TBlioc. doped shows an enhanced su- 
perconducting pairing compared to the parent compound 
model TBlundopod- This model also demonstrates the 
appearance of the (7r,0,0) peak in the non-interacting 
static spin susceptibility. These two features of the 
doped model prove the suggested important role of local 
crystal and electronic structure inhomogeneities due 
to doping for the local superconducting properties of 
Bi2 Sr2 CaCu2 Os+i . 



V. CONCLUSIONS 

A single-band TB model parametrizing the Cu 
3(^2.2 _j^2 bands for the high-temperature superconduc- 
tor Bi2Sr2CaCu208+5 was derived from DFT electronic 
structure calculations. In particular we analyzed the 
changes in the TB model parameters induced by the 
dopant oxygen atom. We found that an accurate quan- 
titative analysis of the dopant-induced changes of the 
electronic structure of Bi2Sr2CaCu2 08+5 requires high- 
quality TB parametrizations. The required quality of 
the modelling was achieved by including effective far- 
neighbor interactions that increased the number of TB 
parameters in the parent compound model up to 14. Two 
possible models for the parent compound were proposed 
and compared, one based on the nearest-neighbors in- 
teractions and the second, based on presumably physi- 
cally justified hybridization paths with an emphasis on 
intra-layer ones. These two parent compound TB mod- 
els were used for constructing the TB models for the 
doped supercell. The nontrivial problem of mapping 
the doped supercell bandstructure was treated by ap- 
proximating the doped supercell TB Hamiltonian either 
by a homogeneous one, with averaged parameters, or 
by one where only certain selected parameters were ad- 
justed to map the DFT bands. The more promising latter 
approach gave results consistent with physical intuition 
when applied to the nearest-neighbors parent compound 
TB model. 

The static spin susceptibility calculated with the doped 
supercell TBlioc. doped model possesses qualitatively new 
features, namely, a pronounced peak at the X = (tt, 0, 0) 
point in the folded Brillouin zone, compared to the sus- 
ceptibility of the parent compound. This change in the 
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susceptibility was shown to lead to a significant enhance- 
ment of d-wave pairing in the same Bi2Sr2CaCu208+i5 
model. Depending on the value of the interaction param- 
eter U chosen, a modulation of the coupling constant A of 
order 30% required to explain the STS phenomenologjPl 
is easy to obtain. While this is not a strictly local calcula- 
tion, it is a strong indication that - within weak coupling 
theory - the local electronic structure change caused by 
the interstitial O dopant in Bi2Sr2CaCu208+5 in fact en- 
hances the pairing locally. A full local calculation of 
the inhomogeneous spin fluctuation pairing interaction 
is clearly desirable. 

Other extensions of the current calculation would be 
useful. Since it is known that the three-band model 
describes the magnetic properties of cuprates better, it 
would be worthwhile to try to extend the approach pre- 
sented in this paper by performing a parametrization of 
the three-band TB Hamiltonian and studying the effects 
of doping on this model's parameters. In this case, one 
expects, on the one hand, a growth of the number of 
model parameters due to new degrees of freedom. On the 
other hand, the number of effective far-neighbor interac- 
tions should decrease as many effects due to hybridization 
between Cu and O orbitals are now incorporated in the 
model itself. One could also consider other multi-hand 
TB models by taking more Cu and O orbitals, besides the 
Cu Sdj.2_y2^ O 2px and O 2py orbitals, into account. The 
multi-band models can be very efficiently parametrized 
within the Slater-Koster formalisnP^, moreover, the for- 
malism allows direct calculation of overlap integrals in 
the doped compound as a function of relative atomic po- 
sitions instead of parametrization by fitting. This will be 
a subject of further investigations. 



tion, the slabs are separated by approximately 15 A of 
vacuum in order to exclude any interaction between the 
marginal Bi atoms of adjacent slabs. The exact atomic 
positions inside the supercell slab were determined by 
He et al^, who performed structural optimization cal- 
culations with the Vienna ah initio simulation program 
(VASP) within the local density approximation. 

Calculations for the parent compound were carried out 
with an energy cut-off for the basis set size given by 
i^MT-K'max = 5.50 (i?MT IS the Smallest muffin tin radius 
and if max is the maximal lattice vector considered) . The 
muffin tin radii for the different atoms in the unit cell 
were the following: i?MT(Bi) = 1.88 bohr, i?MT(Sr) = 
2.22 bohr, i?MT(Ca) = 2.17 bohr, i?MT(Cu) = 1.82 bohr, 
and i?MT(0) = 1.61 bohr. We considered a mesh of 
240 fc-points in the irreducible Brillouin zone (IBZ) cor- 
responding to the space group lA/nimm of the parent 
compoimd unit cell. Both -RiviT-fi'max and the number of 
fc-points in the IBZ were tested to be sufficient for ren- 
dering an accurate electronic bandstructure. 

The supercell has the symmetry of a centered mon- 
oclinic unit cell in the space group Cm which includes 
as symmetry operations the identity and a mirror plane 
reflection perpendicular to the atomic layers passing 
through the interstitial oxygen atom. The atomic dis- 
placements caused by the interstitial oxygen are mirror 
symmetrical about this plane, as illustrated in Fig. |4] for 
Cu atoms. For the supercell calculations, we used the 
same i^Mxi^^max and Rmt values as for the parent com- 
pound and 64 fc-points in the IBZ of the supercell. 

Appendix B: Details of the tight-binding models and 
application to undoped Bi2Sr2CaCu2 08 
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Appendix A: Details of the electronic structure 
calculations 



The Bi2Sr2CaCu2 08+5 supercell is constructed in the 
following way. It is rotated with respect to the primitive 
unit cell by 45° around the 2;-axis and extended along 
the new x- and y-axes such that its xy dimensions are 
(2^20 X 2V2a)i?45°, where a is the x dimension of the 
pure Bi2Sr2CaCu20g unit cell. Then only one slab of the 
two Bi2Sr2CaCu208+5 slabs is considered. In the z direc- 



The parent compound contains four Cu atoms per unit 
cell: two Cu atoms per slab. Since the Bi2Sr2CaCu208 
lattice is base-centered and each eigenvalue in the elec- 
tronic structure is doubly degenerate due to the presence 
of two equivalent slabs, there are only two doubly degen- 
erate Cu 3d^2_y2 bands. We therefore construct the TB 
Hamiltonian by considering only one pair of Cu atoms in 
one slab: 



H = 



k 



[d\{k) 4(k) 



EM 



di{k) 
d2(k) 

(Bl) 

where dj(k) [(ii(k)] create [annihilate] an electron with a 
wave-vector k in the Cu 3dx2-y2 orbital of atom i = 1,2 
(Cu atoms 1 and 2 belong to different Cu02 layers). The 
energy dispersion E^yik) is due to interactions between 
Cu atoms within a layer and Ej_ (k) is due to inter-layer 
electron hoppings; both dispersions are given by 



where ti is the hopping integral between Cu atoms that 
are connected by a vector 1. This vector is denoted as 
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TABLE I; TBlundopcd results: Optimized values of the on-site 
energy jj, and the hopping integrals U between 12 Cu nearest 
neighbors in eV. These parameters reproduce the Cu 3d^2_y2 
bands in Fig. [5] (a) and (b). The vector 1 = (n, m, z) is given 
by integers n, m; z can take values of or z = 0.099 as 0.099c 
is the distance between two Cu02 layers. 





toOz 


iioo tiOz iiio 




0.4212 


0.0543 


-0.5196 0.0056 0.1115 


-0.0221 


^200 




^210 ^212 ^220 ^222 


^300 


-0.0859 


0.0117 


-0.0078 -0.0064 0.0025 -0.0103 


-0.0238 



1 = (n, m, z) and corresponds to a vector (na, ma, zc) 
in absolute coordinates where a and c are the unit cell 
parameters, n, m denote integers and z — 0.099, is the 
distance between two Cu02 layers in units of c. 

By optimizing the values of the hopping integrals ti, 
the eigenvalues of the Hamiltonian Eq. (Bl) can be ad- 
justed to describe the Cu Sd^^. 



_y2 bands of the parent 
compound. In order to obtain a description that satisfies 
our accuracy requirements, we had to include into the 
model effective hopping integrals between 12 Cu near- 
est neighbors as listed in Table [ij We denote this tight- 
binding results TBlundopcd- 

The in-plane nearest-neighbor hopping integral iigo = 
—0.5196 eV and second nearest-neighbor hopping inte- 
gral tiio = 0.1115 eV are in agreement with the results 
of previous DFT calculationa^^l and with the analysis 
of photoemission measurements of the Bi2Sr2CaCu208 
electronic structure^^. In Fig. [2] (a), the energy spec- 
trum of the TB Hamiltonian is compared with the DFT 
bands. The agreement is overall good, except for some 
small deviations around the F point. This region is par- 
ticularly difficult to reproduce since here the shape of 
the Cu 3dx2-y2 bands is strongly influenced by the hy- 
bridizations between the Cu Sd^^^-y^ orbital and other 
energetically close Cu d and O p orbitals. 



Appendix C: Details of the TB model for the 
O-doped supercell 

Here, we present the TB model for the O-doped super- 
cell, based upon the TBlundopcd model for the undoped 
single cell, whose parameters we use as initial. As ad- 
justable parameters, we choose three hopping integrals 
of the tioo type, two of the tuo type and three of the 
^200 typ6 that connect the Cu atoms experiencing the 
largest displacement due to the interstitial oxygen and 
its neighbors (the hopping integrals of the iioo, iiio and 
^200 types are represented by, respectively, solid, dashed 
and dash-dotted lines in Fig.|4]). We also allow for differ- 
ent on-site energies fi for the 3dr^2_y2 orbitals of the 8 Cu 
atoms in the Cu02 layer closest to the dopant. Making 



use of the crystal symmetry, the number of /j, parame- 
ters is reduced to six. We then assign a unique n value 

TABLE II: TB2undopcd results: Optimized values of the on- 
site energy fi and hopping integrals ti for the parent com- 
pound where only the closest nearest neighbors are consid- 
ered. Also shown are the TB results obtained by mapping 
the doped supercell bands to a homogeneous Hamiltonian 
(TB2hom. doped)- The meaning of the three subindices is the 
same as in Table HI 







iioo 


tiio 


^200 


toOz 


TB 2 undoped 
TB2[^om. doped 


0.4464 
0.4900 


-0.5174 
-0.5150 


0.1085 
0.1158 


-0.0805 
-0.0800 


0.0818 
0.0700 






^210 


iaoo 


^400 


t21z 


TB 2 undoped 
TB2hom. doped 


-0.0264 
-0.0229 


-0.0073 
-0.0075 


-0.0182 
-0.0177 


-0.0122 
-0.0046 


-0.0044 
-0.0062 




t220 


t:i30 








TB 2 undoped 
TB2hom. doped 


0.0068 
0.0045 


-0.0052 
-0.0015 


0.0049 
0.0012 


-0.0047 
-0.0003 



to the on-site energies of the other 8 Cu atoms since we 
expect that they are less affected by the dopant. These 
seven on-site energies together with the eight hopping in- 
tegrals are varied during the Hamiltonian optimization. 
The optimized values of the hopping integrals are given 
in Fig. |4] in eV and the on-site energies are listed in the 
Figure caption. We denote this model TBlioc. doped- We 
observe the largest on-site energy variation for the most 
displaced Cu (marked with an arrow in Fig. |4]) while Cu 
further away from the dopant are hardly affected. The 
variations in the hopping integrals are also consistent. 
For example, tioo increases when the two Cu atoms get 
closer and slightly decreases when they are pushed apart 
by the interstitial oxygen. The decrease is even stronger 
when the Cu atoms shift with respect to each other par- 
allel to the mirror plane, what appreciably reduces the 
orbital overlap. 

We also considered the same 15 parameters (7 fi's 
and 8 t's) to map the DFT Cu 3dx2_y2 bands with a 
supercell TB Hamiltonian based upon the TB2undopcd 
for the parent compound. Table [H] While a mapping 
to the doped bands is almost as good as the one given 
by the previous model, the resulting model parameters 
assumed seemingly chaotic values not consistent with 
their expected behavior. One faces similar inconsisten- 
cies also when other trial sets of adjustable parameters 
are used. The failure of the TB2undopcd (Table [h]) in de- 
scribing the dopant-induced changes in the bandstructure 
of Bi2Sr2CaCu208+5 indicates that the results provided 
by the approach based on optimizing certain selected 
model parameters depend very strongly on the choice of 
effective far neighbor interactions that are not optimized. 
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